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Abstract 



'^ , Pinned surface and frozen surface approximations are two commonly used approximations for 

the boundary conditions at the exposed surfaces of semiconductor structures. We have studied 



the effect of pinned surface and frozen surface boundary conditions on the spatial distribution of 
compressible and incompressible strips observed in the two dimensional electron gas formed in a 
GaAs/AlGaAs heterostructure under quantum Hall effect regime. We have used semi classical 
Thomas-Fermi method for describing the many body problem along with the Poisson equation for 



Q ■ electrostatics. We observe that the boundary conditions significantly effect the spatial distributions 



of the compressible and incompressible strips. 
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I. INTRODUCTION 

Quantum Hall Effect (QHE) has continuously been a field of attraction for both theoreti- 
cal and experimental researchers since its discovery.- Theoretical works on this subject were 
concentrated on developing a basis to understand the changes in the longitudinal resistance 
and the formation of plateaus in the transverse (Hall) resistance.-"— These works have led to 
the consensus that the plateaus in the transverse resistance and the changes in the longitu- 
dinal resistance can be understood in terms of compressible strips (CS) and incompressible 
strips (IS). Therefore, determination of the spatial distribution of these strips is important 
for an accurate description of QHE. 

In some experimental works these strips were used in place of electron or photon beams to 
construct devices equivalent to well known interferometer structures such as Mach-Zehnder 
or Aharonov-Bohm interferometers.—"— The spatial distribution of CS and IS is very crucial 
for understanding the behavior of these devices.— i^i^ In another experimental study it was 
observed that changing the magnetic field strength causes a hysteresis in the Hall resistance 
of bilayer systems with a density mismatch.— >2^ It is believed that the underlying physics 
in this hysteresis phenomenon can also be understood in terms of the CS and IS.-'^^i'i^'^^ 
The basic principles for determining spatial distribution of the CS and IS were laid by the 
seminal works of Chklovskiis.-"^ 

CS and IS correspond to partially and fully occupied Landau levels which in turn are 
determined by the applied magnetic field. In a complete description, the charge distribu- 
tion within the 2DEG should be determined from a solution of the Schroedinger equation 
and thermodynamic equilibrium. The effective potential term in the Schroedinger equation 
and the local thermodynamic equilibrium of the 2DEG are determined by the electrostatic 
potential. Therefore, the spatial distribution of CS and IS is basically determined by the 
electrostatics of the structure and a correct description of the electrostatics of the system 
is very important. It is well known that the electrostatics of a semiconductor structure is 
described by the Poisson Equation for dielectric media. The solution of the Poisson equation 
for a specific system is completely determined by the charge distribution and the boundary 
conditions of the system. Thus a correct description of the boundary conditions is very 
important. Naturally the surfaces of the semiconductor structure are the boundaries and it 
is clear that boundary condition for the electrostatic potential is the gate voltage where the 



surface is covered with gates. However, for the exposed parts of the surface the boundary 
condition is not clear. There exists two commonly used assumptions for the boundary con- 
ditions in the literature and they are called as Pinned Surface (PS) and Frozen surface (FS) 
boundary conditions.—""— 

For the PS boundary condition it is assumed that potential on exposed surfaces is fixed 
to the Fermi energy Ep of the 2DEG.— "— Under this assumption exposed surface is con- 
sidered as a constant potential surface at zero potential, or in other words exposed surfaces 
are treated as grounded gates. In fact another plausible assumption for the PS boundary 
condition is that the potential on the exposed surface is fixed to the gate voltage —Vg on the 
gates but in this case 2DEG cannot be confined to the region defined by the gates so that 
assuming the gates to be pinned to zero voltage is more acceptable. This assumption is more 
realistic at higher temperatures such as room temperature. Also, we would like to point that 
this boundary condition corresponds to an unrealistic sudden change of the potential at the 
gate edges from Vg to zero which is unrealistic. Because the value of the potential is fixed 
at the boundary PS approximation corresponds to a Dirichlet type boundary condition.— 

FS approximation corresponds to the case where the charges on the surface are not af- 
fected from the gate potentials or the changes in the charge distribution in the rest of the 
system (i.e. the surface charges are frozen) and these surface charges induce a constant 
electric field in the normal direction of the surface. This approximation is expected to be 
more realistic at low temperatures. The temperatures required to observe QHE is around 
IK so that the FS approximation seems more appropriate for an analysis of QHE. In the 
FS approximation the normal component of the electric field, i.e. the gradient of the elec- 
trostatic potential is fixed on the boundary and this type of boundary condition is called 
a Neumann boundary condition.— It is well known that numerical treatment of partial dif- 
ferential equations with Dirichlet boundary condition is easier than those with Neumann 
boundary conditions. This fact has led many researchers to use PS boundary conditions 
instead of FS boundary conditions.— "—'22i'2^'2£ 

The effects of these boundary conditions on the charge distribution of 2DEG were com- 
pared for some semiconductor structures in the absence of magnetic field.— "^ These works 
have concentrated on the confinement potential on the 2DEG and formation of quantum 
wires or quantum dots. However there is no study comparing the effects of these approxi- 
mations on the spatial distribution of the CS and IS under QHE regime. In this study, we 



will compare the effects of these boundary conditions on the spatial distribution of the IS 
and CS. In order to do so, we use Thomas- Fermi- Poisson (TFP) approach where quantum 
mechanical part of the problem is described by a semiclassical approach, namely Thomas- 
Fermi approximation and the electrostatic part of the problem is described by the Poisson 
equation. TFP approximation has extensively been used for the analysis of systems under 
QHE regime.-"— 

II. THEORY AND THE MODEL SYSTEM 

We consider the system shown in Fig. [H The system was supposed to be infinite and 
invariant in y direction, so that the problem was effectively reduced to a two dimensional 
problem. Also we have assumed that the structure is periodic in the x direction and we have 
used periodic boundary conditions in this direction. The z direction in Fig. [1] corresponds 
to the growth direction. We have assumed that the structure has a back gate at V, i.e. 
the substrate is grounded. So that the boundary condition at the bottom is Vz=5^m = 0. 
The boundary condition at the top surface is V^=o = ~Vg at the gates. Other parts of the 
top surface of the structure are exposed surface. We have performed separate calculations 
using both the FS and PS boundary conditions for this region. 

Following Davies et al.— "— , we have taken the electrostatic potential to be zero on the 
exposed surfaces under the PS approximation. The top surface separates the space into two 
domains and Poisson equation can be solved separately in each domain. However in this 
case we do not need the solution of the Poisson equation in the domain filled by air and 
solving Poisson equation for only the domain filled by the semiconductor structure in Fig. [T] 
is sufficient. 

In the FS case we have assumed that the charges cannot move within the semiconductor 
structure so that the potential on any surface is not pinned to any constant value. In this 
case we treat the exposed surface as a free boundary and we assume that the gradient 
of the electrostatic potential is zero at some point sufficiently far from the surface of the 
structure. There are different layers within the semiconductor structure each one having 
different dielectric constants. However the differences between these dielectric constants are 
small and we have assumed that the dielectric constant is the same all over the semiconductor 
structure. However, considering the significant change in the dielectric constant on the 
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FIG. 1. A sketch of the model heterostructure used in the calculations. 



exposed surface we impose continuity conditions on both the electrostatic potential and the 
dielectric displacement. 

A sketch of the structure considered in this work is given in Fig. [1] The depth of the 
semiconductor structure is taken as 5 fim along z direction. For the FS calculations, the 
healing distance in the air is taken as 3 fim. The extent of the semiconductor structure 
along X direction is 2 fim. Each gate is assumed to be 150 nm wide, but due to the 
periodic boundary conditions along x-direction the effective width of each gate is 300 nm. 
We have assumed that (5-doping is used in fabrication of the structure. The depth of the 
donors from the surface is taken as 50 nm while 2DEG is at 90 nm from the surface. The 
surface density of ionized donors and average electron density were assumed to be equal 
nd = fie = 3.0 X lO^^m"^. 

Encouraged by the slowly changing background potential, we have used the Thomas- 



Fermi (TF) approximation for the description of the quantum mechanical part of the prob- 
lem. TF approximation presents a semi classical approach to obtain electron density in 
many-electron systems; it was used efficiently in many studies relevant to QHE because it 
provides simplicity in the calculations.-"- Within TF approximation the electron density is 
given as 

n,i{r) = fdED{E)f[{E + V{f) - ^)/A;bT] (1) 

where D{E) is density of states, and V{r) is the electrostatic potential, fi is the chemical po- 
tential, kB is Boltzmann constant and T is the temperature while f{a) = {l+exp{a))^^ is the 
well known Fermi distribution function. In the absence of magnetic field the density of states 
for 2DEG is constant and given as D{E) = m/Tifi^. On the other hand, under the influence 
of a magnetic field of strength B the density of states becomes D{E) = ^^ J2'^=o ^(-^ ~ -^") 
here f^s = 2 is spin degeneracy, i = ^/fijeB is magnetic length. En = (n -|- 1/2) ^c is energy 
of nth Landau level and Uc = eB/m is cyclotron frequency. Therefore the electron density 
is given by 
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in the absence and presence of a magnetic field, respectively. 

The chemical potential yU appearing in Eqs.(2), is determined by the condition of conser- 
vation of the total number of electrons, that is 

No= fd^fUelif) (3) 

here A'^o is total number of electrons. On the other hand the electrostatic potential must 
satisfy Poisson equation so that 

V. (eVV(f)) = -p(f) (4) 

where plf) is total charge density and e is the dielectric constant of the material. Total 
charge density can be written as p{r) = Pdir) — enei{Vr) where Pd{r) is the fixed charge 
density due to donors and impurities except 2DEG. Equations (2-4) form a closed set of 
coupled equations. Electrostatic potential, chemical potential and electron distribution in 
the 2DEG should be obtained from a simultaneous solution of these equations. However, 



during the simultaneous solutions of these equations one may encounter many divergence 
and instability problems. In order to overcome these problems we have used a slightly 
modified version of an iterative algorithm previously used by Gerhardts and co-workers.-~- 
This algorithm can be summarized as follows: 

First, in order to create an initial value for the iterative steps, Eqs.(2-4) are solved self- 
consistently and electron distribution is obtained in the absence of magnetic field at absolute 
zero temperature. Then using this first solution as an initial value at a sufficiently high 
temperature T in the absence of magnetic field Eqs.(2-4) are solved again self-consistently 
and chemical potential /x, electrostatic potential V{f) and electron density nei{r) are ob- 
tained. We use this finite temperature solution as an initial value to obtain a solution for 
the Eqs.(2-4) for the same temperature but now in the presence of the magnetic field. Af- 
terwards the temperature is gradually reduced until we reach the target temperature. For 
the calculations presented in this study, we have chosen the initial high temperature and 
the target temperature as 60 K and 1.4 K, respectively. While decreasing the temperature, 
the result obtained in the previous step is used as the initial value at each step. For the 
iterative solution, which forms the basic step of this algorithm, we used a generalized form 
of Newton-Raphson Method details of which is described below. 

Let fi* and V*{f) be the chemical potential and the electrostatic potential corresponding 
to the exact simultaneous solution of Eqs. (2-4). On the other hand we use the potential 
function V{f) obtained in the previous step as a first estimate. If we indicate the difference 
between the exact and estimated potentials as SV{r), then we can write the exact potential 
as 

V*{f) = V{r) + 6V{r). (5) 

The exact potential satisfies Poisson equation, i.e. 

VV* = VV + V^6V = --{pd - erieiiV + 6V)). (6) 

If we expand riei in a Taylor series about V and keep only the first order terms 

VV + V'6V = -i(prf - en,i{V + 6V)) + ^^SV (7) 

is obtained. According to Eqs. (2), we can write the relation between derivatives of electron 
density riei with respect to electrostatic and chemical potentials as 

duei duei 



dV dn 
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(8) 



So Eq. (7) can be written as 

VV + V^5V = -^(p, - en,i{V + 5V)) - -^^SV (9) 

or 

Here if we define the functions p = ^^^ and / = —V'^V — \{pd — en^iiV + 5V)) then Eq. 
(10) reduces to 

V'^SV + p5V = f. (11) 

Equation (11) is a second order differential equation for the correction term 5V however it 
involves a linear term unlike the Poisson equation. By adding 5V , which is obtained from 
the solution of Eq. (11), to the previous potential estimate we obtain a better estimate for 
the electrostatic potential function. If we denote the potential estimates obtained in two 
successive steps as V^"+^ and V^ then one step of the Newton-Raphson algorithm can be 
expressed as 

Then the chemical potential and electron density which correspond to this new potential 
are obtained by using Eqs. (2) and Eq. (3). This procedure is repeated until the changes 
between two successive steps in the potential, electron density and chemical potential are 
small enough. 

III. RESULTS AND DISCUSSION 

It is already pointed out that the boundary conditions are very important in obtaining 
the solution of Poisson equation. These boundary conditions refiect the structure and the 
behavior of the surface and have important effects on the behavior of 2DEG. In fact the 
purpose of this study is to investigate how different boundary conditions affect the behavior 
of 2DEG under intense magnetic field. Therefore we have performed calculations for phys- 
ically plausible two different boundary conditions for the upper (exposed) surface of the 
structure. These boundary conditions correspond to two different surface behaviors, namely 
PS and FS approximations. We have used the model structure shown in Fig. [1] in all our 
calculations. We have chosen the gate voltages so that the potential difference between the 
gates and the center of 2DEG is —0.35 V . 

8 



We present our results in Fig. |2] for electron density distribution and variation of the 
screened potential under both FS and PS boundary conditions for a magnetic field intensity 
of 6.0 T. The solid(black) lines correspond to the PS and the dashed(red) lines correspond to 
the FS boundary conditions. In order to compare the potentials we have given the potential 
difference from the edge of the structure in the plane of 2DEG. Under the PS boundary 
conditions, the potential at the top surface of the structure is forced to an unphysical jump 
at gate edges, even though this sharp change smooth as we move away from the surface 
its effect still can be seen in the plane of 2DEG. In contrast to the PS approximation, the 
potential in the FS approximation is not forced to a jump at the exposed surface, and thus 
the potential profile on the plane of 2DEG also exhibits a more smooth behavior. The 
sharp change in the potential results in a wider depletion region (i.e. a region without any 
electrons) at the edges and rises more quickly than those obtained under FS approximation. 
The smoother change in the potential under FS approximation results in a slower change in 
the electron density and under this approximation electrons penetrate more toward the edge 
of the structure so that the depletion region is narrower. While the electron density changes 
more rapidly near the edges under PS approximation, it quickly saturates and behaves like 
constant near the center. This quick saturation should be expected because the distance 
between the exposed surface and 2DEG is very small and the potential is assumed to be 
constant at the exposed surface. This behavior of the electron density changes completely 
for FS approximation where the electron density and the electrostatic potential changes 
more slowly. 

Because the structure is symmetrical with respect to the center and the changes in the 
potential and the electron density occurs near to the edges we plot the rest of our results for 
a smaller range at the left edge of the structure. In Fig. [3] we present the electron densities 
obtained for various magnetic field strengths. The potential profiles for the same magnetic 
field strengths are given in Fig. |H It can be seen from Fig. |3]that IS move toward edge of the 
structure and get narrower with decreasing magnetic field. Hall plateaus in the system occur 
if there is at least one IS in the systen>22i^ and when the width of the IS is less then cyclotron 
radius of an electron, the system leaves the plateau.- When one uses the PS approximation 
the first appearance of IS occurs at lower magnetic field strengths. In Fig. |3l there is no 
IS observed for B = 7.75 T and IS appears for B = 7.5 T under PS approximation while 
under FS approximation IS is observed for both B = 7.75 T and B = 7.5 T. Also the sharp 
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FIG. 2. Variation of the electron density and screened potential in the 2DEG for B = 6.0 T.The 
solid(black) lines correspond to the PS and thB^ashed(red) lines correspond to the FS boundary 
conditions. 
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FIG. 3. Spatial electron density distributions in terms of filling factor for several magnetic fields. The 
solid(black) lines correspond to the PS and the dashed(red) lines correspond to the FS boundary 
conditions. 

change of the potential under PS approximation results in even narrower IS. Thus under 
PS approximation one predicts the system to leave plateau region at a higher magnetic field 
strength. As a result, using PS approximation may lead to wrong predictions of narrower 
plateau regions. 
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FIG. 4. Spatial distribution of the screened potential for several magnetic fields. The solid(black) 
lines correspond to the PS and the dashed(red) lines correspond to the FS boundary conditions. 

The total electrostatic potential is called as screened potential when it is obtained self 
consistently because it also contains the screening effect of the 2DEG. Screened potential 
profiles obtained for various magnetic field strengths are given in Fig. |H Because in the 
CS region 2DEG behaves like a metal the potential is perfectly screened and behaves like 
constant. Whereas in the IS regions the potential cannot be screened and changes until it 
reaches to the next Landau level. So the screened potential is constant while the electron 
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density changes and electron density stays constant while the potential changes. However, 
the effect of the two approximations for the boundary conditions at the exposed surface on 
screened potential profiles is similar to the electron density. 

It is generally accepted that FS boundary conditions reflect the real physical conditions 
better for semiconductor structures at low temperatures.—'^'^ However, one has to consider 
the region beyond the exposed surface in some way under FS boundary conditions. When 
working with FS boundary conditions a complete description of the problem requires the 
solution of Poisson equation for each electron distribution in the 2DEG. On the other hand, 
working with PS boundary conditions is much easier. The exposed surface under PS bound- 
ary conditions is treated as a grounded gate and the mirror charge convention may be used 
for a complete description of the problem. Then solving the Poisson equation for the gates 
and the donor distribution only once will be sufficient. In conclusion, we observe that using 
PS approximation for exposed surface boundary conditions leads one to obtain erroneous 
results for the spatial distribution of IS and CS under QHE regime. However most of the 
studies in the literature refrain from using FS boundary conditions because it is easier to 
work with PS.^-iS'2a^'2^ 
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